UHE tau neutrino flux regeneration while skimming the Earth 
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The detection of Earth-skimming tau neutrinos has turned into a very promising strategy for the 
observation of ultrahigh-energy cosmic neutrinos. The sensitivity of this channel crucially depends 
on the parameters of the propagation of the tau neutrinos through the terrestrial crust, which 
governs the flux of emerging tau leptons that can be detected. One of the characteristics of this 
propagation is the possibility of regeneration through multiple v T <-> r conversions, which are often 
neglected in the standard picture. In this paper, we solve the transport equations governing the 
v T propagation and compare the flux of emerging tau leptons obtained allowing regeneration or 
not. We discuss the validity of the approximation of neglecting the v T regeneration using different 
scenarios for the neutrino-nucleon cross sections and the tau energy losses. 



I. INTRODUCTION 

With the advent of a new generation of large-scale 
detectors of cosmic radiation, the observation of high- 
energy cosmic neutrinos produced in distant astrophys- 
ical sites or possibly by other, more exotic, mechanisms 
has become one of the major challenges of astropar- 
ticle physics. In both astrophysical and exotic mod- 
els, substantial fluxes of electron and muon neutrinos 
are expected from the disintegration of charged pions 
(and kaons) produced in the interaction of accelerated 
particles with ambient matter and radiation, either at 
the source location or between the source and the ob- 
server. Given the large distances traveled by these cos- 
mic neutrinos, approximately equal fluxes in v e ,v^ and 
v T are expected on Earth as a result of flavor mixing 
and oscillations [l], Q • Important efforts are ongoing to 
build dedicated neutrino telescopes both in the South- 
ern [!, 0] and Northern [E [E @, [1] hemisphere, open- 
ing a window on the neutrino sky in the energy range 
10~ 6 EeV < E v < 10 _1 EeV. At even higher energies, 
other promising experiments are developing the detection 
of coherent radio emission produced by neutrino-induced 
showers in matter 
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Interestingly, recent studies have shown that the new 
generation of ultrahigh-energy (UHE) cosmic ray detec- 
tors such as the Pierre Auger Observatory [16[ and the 
HiRes Fly's Eye detector [17| have a comparable detec- 
tion potential for UHE neutrinos in the range of ener- 
gies 10 _1 EeV < E v < 10 2 EeV, where neutrinos are 
expected to be produced in the interaction of UHE cos- 
mic rays with the cosmic microwave background [181 ] . It 
is long known that downward-going showers induced by 
neutrinos that penetrate deep in the atmosphere can in 
principle be identified at large zenith angles (9 > 75°) 
where no background is expected from hadronic pri- 
maries [l9|,[2(|. It was also pointed out more recently that 
the presence of v T in the cosmic neutrino flux provides 
another promising channel of detection for air shower de- 
tectors [2l|, [22| . Upward-going UHE tau neutrinos that 



graze the Earth just below the horizon (often referred to 
as "Earth-skimming neutrinos" ) are indeed likely to in- 
teract in the crust and produce a tau lepton which may 
emerge and initiate an observable air shower, provided it 
does not decay too far from the detector. 

The sensitivity to such UHE Earth-skimming neutri- 
nos crucially depends on the conditions of the v T prop- 
agation through the terrestrial crust, and on the correct 
estimation of the flux of t leptons that emerge from the 
Earth. This propagation problem has been widely dis- 
cussed in different contexts and with different approxi- 
mations [H, H HE HE S3, HI S<| . An exhaustive treat- 
ment should account for r and v r neutral-current (NC) 
and charged-current (CC) interactions with nucleons, r 
decay and energy losses. However the full coupled trans- 
port equations admit no analytical solution, and even in 
the case of Monte Carlo calculations, simplifications are 
usually made such as dropping the r weak interactions 
and neglecting multiple regenerations of the v r . 

Such approximations are defendablc in the standard 
case where the characteristic lengths for the r CC inter- 
action (~ 600 km, at 1 EeV) and for the r decay (<~ 50 
km) are larger than that for energy losses (~ 6 km). How- 
ever, they might be challenged in other, more exotic sce- 
narios. The knowledge of the neutrino cross section and 
the tau energy losses in the EeV energy range is indeed 
limited [IE [3l[ , and these could be significantly affected 
at center-of-mass energies beyond the TeV by the onset 
of new physics beyond the standard model. Several stud- 
ies have even suggested that the comparison of the flux 
rates between down-going and Earth-skimming neutrinos 
could actually help constraining the neutrino properties 
at ultrahigh energies, where no direct measurements ex- 
ist [Hili. 

All these considerations pinpoint the necessity of an 
accurate determination of the flux of emerging r leptons. 
In this context, it is important to correctly describe all 
the contributions to the r flux and to assess the impact 
of simplifications in the description of the propagation 
problem. 

The present paper therefore focuses on understanding 
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FIG. 1: Geometry of the transport problem. 



the effects of the v r and r regeneration while skimming 
the Earth on the flux of emerging r leptons in different 
scenarios for the neutrino-nucleon cross section and for 
the r energy losses. In Sec. HH we present the transport 
equations for the v and r propagation, and the scenarios 
to be studied. In Sec. IIII1 we compare the flux of emerg- 
ing taus for different conditions of propagation. Finally, 
we present our conclusions in Sec. IIV1 



II. 



v T PROPAGATION THROUGH THE EARTH: 
THE GENERAL PICTURE 



The geometry of the propagation problem is described 
in Fig[TJ where an example of regeneration chain through 
multiple CC interactions and r decays is sketched. Given 
a beam of parallel neutrinos incident on the Earth at a 
given angle a, the problem becomes unidimcnsional and 
the flux of r leptons that emerge only depends on the 
path length traveled across the rock. In the hypothesis 
of a spherical Earth with a crust of constant density, this 
length is directly related to a, which is also the angle of 
the emerging tau. 



A. Transport equations 

We reproduce here a classical formulation in terms of 
the transport equations that describe the evolution of the 
v T and t fluxes, <1>„ T and (f> T , along their path through 
the Earth, accounting for all possible production and ab- 
sorption processes taking place within an infinitesimal 
da;: 
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Here A^ c and A^ 7 are the mean free paths corre- 
sponding respectively to NC and CC interactions of the 
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FIG. 2: Flux of emerging r's (scaled by E 2 ) obtained by 
solving the transport Eqs. fT} and for incident tau neu- 
trinos with an angle a up to 15° and injection flux dN/dE — 
4.6 x I0 7 E~ 2 EeV -1 sr -1 . The points correspond to the 
Monte Carlo solution, while the line is the solution obtained 
with the iterative method. 



incident v T , while a^ c corresponds to the r CC interac- 
tion, which regenerates a v T . T T is the tau lepton lifetime. 

The transport equation for the r can be inferred in a 
similar way. In this case, the small inelasticity of the r 
radiative interactions allows us to use the continuous en- 
ergy loss approximation [34]. Defining ■j(E) = — dE/dx, 
the equation reads: 

d$ T (E,x) ® T {E,x) <f> T (E,x) <*> T (E,x) 



dx 



+ -^h(E)<S> T (E,x)) 



+pATa 



\dec( E ) \?C(E) X? C {E) 
^? c (y,^- y ) 



1-y ' \l-y' J dy 
1-y Ur \l-y' J dy 



(2) 



where X^f and \^ are now the mean free paths associ- 
ated with the r NC and CC interactions, and A^ ec is the 
decay length corresponding to T T . 

We are thus left with a system of two coupled, integro- 
differential equations, which cannot be solved analyti- 
cally without making simplifying assumptions on the na- 
ture and the respective importance of the interactions 
that both the v T and the r can undergo. To preserve 
the generality of the solution we use a Monte Carlo sim- 
ulation that includes all the processes listed above and 
follows the incident particle all the way through the rock, 
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deciding at each step on its fate according to the distri- 
butions encoded in Eq[T] and EqO 

An iterative method that also allows us to avoid sim- 
plifications was proposed in [35|, and further discussed 
in I27l| a nd in [36|. Our equations differ from those 
in [23,[3y| in that we also include the two terms describing 
r NC interactions in the transport equation for $ T . Tak- 
ing into account this generalization, we have applied the 
iterative method (see the Appendix for details) to cross- 
check the results from the Monte Carlo simulation. In 
Fig. we show the good agreement between both calcu- 
lations for the flux of emerging t leptons given an incident 
tau neutrino flux dN/dE = 4.6-10 7 E~ 2 EeV" 1 sr" 1 . 



B. Cross Sections 

In the framework of the standard model, the neutrino- 
nucleon (i/N) CC and NC cross-sections, and crf? c , 
describe deep-inelastic scattering processes. They are ex- 
pressed in terms of the structure functions of the nucleon, 
which in turn depend on the individual parton distri- 
bution functions (PDFs). The PDFs are obtained from 
measurements at accelerators, in determinate ranges of 
Bjorken-x and momentum transfer Q 2 . However, the 
range of parameters probed by the UHE neutrinos, x ~ 
1CT 5 (1 EeV/E v ) and Q 2 up to - 1CT 5 EeV, is outside 
the measured domain and therefore extrapolations are 
needed. 

We choose as a benchmark for the standard i>N cross 
section a recent parametrization presented in [30j , which 
fits to updated HERA data: 



a^(E) = 2Aa^(E) 



6.04 ■ 
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As pointed out in the same paper, another, more spec- 
ulative approach based on the color glass condensate for- 
malism [33] has also been put forward recently to account 
for saturation at very low x. We will also consider this 
case and use the following approximate parametrization 
deduced from Fig. 1 of [30]: 



7 low( E ) 



= 2.4 a. 



= 3. 



low [e) 



t? v 0.170-0. 037-log 10 • 

EeV / 



10 3 pb (4) 



as an example of low vN cross section. 

On the other hand, plenty of models using new physics 
predict an enhancement of the neutrino cross section [38| . 
As an example of high uN cross section, we will use 
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FIG. 3: Flux of emerging taus corresponding to incident tau 
neutrinos with an angle a up to 3° or 15° and dN/dE = 
4.6 x 10 7 E~ 2 EeV" 1 sr 1 using the standard combination of 
parametrizations a st d ® Pstd as defined in Sees. Ill Bl and HTCl 
The regeneration is neglected in the dashed histograms and 
included in the solid ones. 



C. Tau energy losses 

The t energy losses are usually parametrized using the 
following formula: 



dE T 
~dX 



= -a- (3{E T ) E T . 



(0) 



where a is practically constant and accounts for ioniza- 
tion processes, while 0{E T ) includes the radiative contri- 
butions from bremsstrahlung (b), pair production (pp) 
and photonuclear (pn) interactions of the r. 

For the purpose of our analysis, we adopt the same 
strategy as in the previous section and extract from exist- 
ing computations [U Hi| a panel of ad- hoc parametriza- 
tions of (3(E T ) in the relevant range of energy, represen- 
tative of the values found in the literature. The set of 
parametrizations under study is the following (valid in 
the range 10" 1 EeV < E v < 10 3 EeV): 



Pstd(E T 



1.2 + 0.16 x ln( 



E T 



10 EeV 



Phigh(E r ) = 1.36 C 



p s 0.35 

EeV/ 



1 cm 2 



01 OW (E T ) = (^0.3 + 0.06 xlog 10 (^)J lO^g-W 

respectively for a standard, a high and a low value of the 
energy loss parameter (3. 



III. THE FLUX OF EMERGING TAUS: 
RESULTS AND COMPARISONS 



aZ c gh (E) = 2.4 a^ c gh {E) = 3 • ag c d (E) (5) 
without assuming any particular model. 



The Earth-skimming technique has been recently used 
to search for v T with the HiRes telescopes 17[ and the 
Pierre Auger Observatory [Jfj. While both detectors 
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FIG. 4: Flux of emerging taus produced by incident <& Vt (E v ) 
= 4.6- 10 8 sr -1 with energy E v = 0.3, 3 and 30 EeV (from 
top to bottom) and an angle up to 3°. The standard combi- 
nation of parametrizations a a td ® Pstd as defined in Sees. Ill Bl 
and lll Cl is used. The solid histogram takes into account all re- 
generation channels. The dotted histogram only accounts for 
the regeneration through the t NC interaction. The dashed 
histogram neglects the regeneration. 



have a similar energy threshold (~ 0.1 EeV), they are 
sensitive to a different range in the emerging angle a of 
the t leptons. The latter detects the particles from the 
air shower initiated by the r decay that reach the ground 
and, hence, the detector will only be triggered by r lep- 
tons slightly up-going (a < 3°). The former is designed 
to record the faint ultraviolet light emitted by nitrogen 
molecules that are excited as the shower traverses the at- 



FIG. 5: Flux of emerging taus produced by incident <& Vt (E v ) 
= 4.6T0 8 sr" 1 with energy E v = 0.3, 3 and 30 EeV (from 
top to bottom) and an angle up to 3°. The combination 
<?high 8> Pstd as defined in Sees. Ill Bl and ITTCl is used (with the 
same histogram code as in Fig. [4]| . 



mosphere; it is therefore sensitive to t leptons emerging 
with larger angles. 

In figure [3l we present the results obtained with the 
Monte Carlo approach for the flux of emerging t leptons 
with an angle a up to 3° and 15° produced by an incident 
tau neutrino flux dN/dE = 4.6 10 7 E~ 2 EeV" 1 sr" 1 . 
All cross sections and energy loss parameters are here 
set to standard values, as defined in Sees. Ill Bl and III CI 
The effect of the regeneration is negligible if only almost 
horizontal emerging r leptons (a < 3°) are taken into 
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FIG. 6: Flux of emerging taus produced by incident <& Vt (E v ) 
= 4.6- 10 8 sr -1 with energy = 0.3, 3 and 30 EeV (from 
top to bottom) and an angle up to 3°. The combination 
o"std <8> Plow as defined in Sees. ITTBh nd [TT"Cl is used (with the 
same histogram code as in Fig. [4}. 



FIG. 7: Flux of emerging taus produced by incident <& Vt (E v ) 
= 4.6-10 8 sr" 1 with energy E v = 0.3, 3 and 30 EeV (from 
top to bottom) and an angle up to 3°. The combination 
&high (Si Plow as defined in Sees. Ill Bl and HTCl is used (with the 
same histogram code as in Fig. [4| . 



account. On the other hand, if the detector is sensitive 
to larger angles and energies down to a fraction of EeV, 
neglecting the regeneration becomes important already 
within this standard picture (30% less r leptons at 0.3 
EeV for an incident neutrino flux dN/dE oc E~ 2 ). 

The approximation of neglecting the regeneration 
seems thus to be safe in the standard picture for detectors 
only sensitive to almost horizontal showers. However, we 
will now show that this approximation is questionable 
even for them as soon as one considers non-standard sce- 



narios. We have investigated the effect of the regenera- 
tion for those detectors using the set of parametrizations 
defined in Sees. Ill Bl and III CI These expressions do not 
constitute an exhaustive list of all the alternative and 
exotic models that exist in the literature, but rather a 
selection of a few examples that allow us to point out 
the relevance of this process. The combinations of cross 
sections and tau energy losses should be chosen on basis 
of coherent PDFs but this coherency is hard to assert 
even inside the frame of the standard model. Therefore, 
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we chose to investigate the regeneration effects for the 
3® 3 possible combinations, which is sufficient to provide 
a qualitative answer to the question raised in this paper. 

As a starting point, one should however notice that the 
regeneration requires the r to be converted back into a 
v T before losing too much energy, and then this v T to un- 
dergo another CC weak interaction to produce a r again. 
Hence, one expects that a higher energy loss or a lower 
cross section will reduce the effect of the regeneration, 
while, on the other hand, a higher cross section or a lower 
energy loss will enhance it. This is indeed what we find. 
In the following, we focus on the particular combinations 
for which the regeneration is not negligible anymore, and 
compare them to the standard case, essentially on the 
basis of monoenergetic incident beams of i/ T 's. This is 
because the E~ 2 power-law spectrum usually assumed 
for cosmic neutrinos, although adequate to look at the 
global flux of emerging r leptons, is completely domi- 
nated by the lowest energies, while the behavior at the 
highest energies may be important for harder fluxes of 
incident neutrinos. 

In Fig. HE we show the flux of r leptons produced by 
incident neutrinos at given energies (0.3, 3 and 30 EeV) 
and with an angle a up to 3°, for the following combina- 
tions of parametrizations: a st( i <8> (3 s td (the same as used 
in Fig. |3J), (Thigh ® Pstd, °std ® Plow and <r hlgh ® f3i ow . For 
a standard choice of parametrizations, one can see from 
Fig.H]that the impact of regeneration is completely negli- 
gible at low neutrino energies. As the v r energy increases, 
however, the lowest-energy part of the r spectrum starts 
to be significantly underestimated (15% less taus for 30 
EeV incident neutrinos). This effect is washed out in 
the case of an incident neutrino flux dN/dE oc E~ 2 , 
for which the contribution of the highest-energy part of 
the neutrino spectrum is negligible; but it could affect 
the spectrum of emerging r's if the neutrino flux were 
harder. For the other cases, the lowest-energy range of 
the r spectrum becomes underestimated already for in- 
cident v T with energy of 0.3 EeV. For the combination 
Thigh <8> l3i ow , one loses about 70% of the r leptons emerg- 
ing from the Earth if the regeneration is neglected, even 
if assuming an incident flux of neutrinos dN/dE oc E~ 2 
(Fig.©. 

The neutrino flux can be regenerated through the r 
decay and the r CC weak interaction. We present our 
results in a form that allows us to disentangle the two 
effects: the dotted histograms in Fig. I4I5I6I7I account for 
regeneration through the weak interaction of the r only, 
while the solid ones include both effects. As can be ap- 
preciated in the figures, the CC contribution is negligible 
at all energies in most scenarios. It only becomes impor- 
tant for 2/ T 's at the highest energies (~ 30 EeV), in the 
energy range between 0.3 and 3.0 EeV of the emerging 
r spectrum, if a high cross section is assumed (Figs. [3] 
and [7]). There, the effect is at the level of that from the 
regeneration through the r decay. 




FIG. 8: Flux of emerging taus corresponding to incident tau 
neutrinos with an angle a up to 3° and dN/dE = 4.6 x 10 7 
E- 2 EeV" 1 sr 1 using the combination of parameterizations 
(Thigh® Plow as defined in Sees. Ill Bl and llTCl The regeneration 
is neglected in the dashed histogram and included in the solid 
one. 



IV. CONCLUSION 

In this paper, we have studied the mechanism of re- 
generation of the v T flux while crossing the Earth, to 
investigate its effect on the flux of emerging r leptons. 
Assuming a detector with an energy threshold of 0.1 EeV 
and sensitive only to taus with an emerging angle below 
3°, the effect is negligible for a flux of incident neutrinos 
dN/dE oc E~ 2 with the standard values of cross sections 
and energy losses. But we have shown that this is not 
valid for other assumptions on the detector performance, 
neither for less standard values of the weak cross section 
or tau energy losses. On the one hand, neglecting the 
regeneration leads to about 30% underestimation on the 
integrated flux for a detector sensitive up to 15°. On the 
other hand, the underestimation of the r's emerging from 
the Earth with an energy greater than 0.1 EeV and an 
emerging angle below 3° may reach 70% for the different 
values of a Vr and /3 T studied in this paper. Moreover, 
the error made in the computation of the flux of emerg- 
ing taus can be even larger for a harder flux of incident 
neutrinos since the contribution of the regeneration effect 
increases with the energy of the incident neutrino. 

The simplification of neglecting the regeneration is 
thus only safe for particular values of the physical proper- 
ties playing a role on the propagation and specific detec- 
tors. It may lead to a significant underestimation of the 
flux of emerging r's when looking at nonstandard values 
of the weak cross section or tau energy losses. There- 
fore, it should be carefully treated and accounted for 
when studying the systematics due to the uncertainties 
on those properties or while using the Earth-skimming 
technique to test for instance higher weak cross-sections. 
Similarly, one should carefully check the effect for the 
characteristics of the actual detector before neglecting 
the regeneration. 
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APPENDIX A: ITERATIVE SOLUTION OF THE 
TRANSPORT EQUATIONS 

We reproduce here a numerical method proposed in 
(27l [35| , generalizing it to take into account the NC weak 
interaction of the r lepton. We start from a solution of 
order and iterate it up to get a solution that satisfies 
the relation F($(E, x),d x $(E, x)) = as given by Eq.Q] 
or [2] For numerical reasons, it is better to work with 
numbers between and 1 so that we write the equations 
in terms of the deviation to the solution of order 0, $° , 
rather than in terms of $> Vt itself. We therefore assume 
that the solution is of the form 



$„ T (£,z) = $° T (£)exp 



K T {E) 



(1-Z Vt (E,x)) 



where X^(E)- 1 = A£ c (E)- 1 + A^ c '(E)- 1 accounts for 
the possible sources of depletion of <&„ T at an energy E, 
and Z Vt (E, x) is a function to be determined, which in- 
cludes all contributions to the regeneration of <I>„ T at the 
same energy. 

Z Vt (E, x) can be obtained iteratively by forcing it to 
satisfy the transport equations, viewed as a recurrence 
relation with the initial condition Z® (E, x) = 0: 



z:^(E,x) = -. 



dZ^(E,x) 



dx 
dy 



+pAfA $uE, X )j i-irni-y" 



E 



-pM A 



dy 



®Z T (E,x) J 1-y T \l-y 



dy 
rCCr„. E 



E ^da^iy,^) 



dy 



+ cS>»(E,z)J l-y^\l-y' X J dy 



As the two equations are coupled, it is necessary to 
know $ T at a given order to proceed with the iteration at 
the next order. Following the same philosophy as above, 
we assume that the solution of Eq[2]for $ T is of the form: 



§ T {E,x)= duS(E(u-x),u)K(E,u,x) 
Jo 

where the following notations have been introduced: 



S(E,x) = P Ma 



K(E 1 u,x) = exp 



l-y ^\l-y 



a*g c (y,^- y ) 



dy 



k{E,u,x) ■ (1 - Z T (E, x) 



k(E, u, x) 



+ Z T {E(u-x),u)) 

~d~/(E(v-x)) 



dv 



dE 



Xr(E(v-x)) 



'AiidKiE)- 1 = Xf^E)- 1 + \C C [E)- 1 + \* c (E)- 1 . The 
equation for ZV: then reads: 



x) 



x / duS n+1 (E{u- x),u)K(E,u,x) 



Z?(E(u - x), u) - k(E, u, x) 



d x Z?(E,x) 



(d El {E)-K\E)) 



{d El {E)-\r'(E))^ T {E,x) 



again with Z®(E,x) = 0). Practically, the iteration 
can be stopped as soon as n=3. 
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